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C^) ' We have developed an efficient method for quantum number projection from most gen- 

eral HFB type mean-field states, where all the symmetries like axial symmetry, number 
conservation, parity and time-reversal invariance are broken. Applying the method, we 
have microscopically calculated, for the first time, the energy spectra based on the exotic 
tetrahedral deformation in 108,110 Zr. The nice low-lying rotational spectra, which have all 
characteristic features of the molecular tetrahedral rotor, are obtained for large tetrahedral 
deformation, 032 > 0.25, while the spectra are of transitional nature between vibrational 
and rotational with rather high excitation energies for 032 ~ 0.1 — 0.2. 
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§1. Introduction 

The progress of radioactive beam facilities in these years provides us a great 
possibility to explore nuclear regions with various combinations of neutron and pro- 
ton numbers, where many interesting new phenomena have been predicted. In the 
present work, we would like to focus, among others, on the exotic shape with high 
rank point group symmetry j-Q i.e., the tetrahedral nuclear states. This "pyramid- 
like" shape leads to the extra stability in shell energy due to the higher symmetry, 
e.g., the appearance of fourfold degenerate orbits. It has been suggested that such 
states appear as low-lying states, or even as a ground state, in some nuclei around the 
"tetrahedral-closed shell" nuclei® see also Refs. E]),H]),[S]),ED- In these earlier studies, 
the mean-field approaches are employed to search for minimum configurations with 
the tetrahedral shape. However, the quantum excitation spectra and the properties 
of electromagnetic transition rates between them are necessary to definitely iden- 
tify the tetrahedral states. The quantum number projection, especially the angular 
momentum projection, from mean-field states is useful for such a purpose. 

Note that the tetrahedral shape is neither axially symmetric nor reflection sym- 
metric; the number conservation is lost if the pairing correlation is included, and the 
time-reversal invariance is broken if the cranking procedure is used to investigate 
the collective rotational bands. Therefore, the most general HFB type mean-field 
states without symmetry restrictions are necessary to describe the tetrahedral nu- 
clei. Recently we have developed an efficient method for quantum number projection 
from such most general HFB type states.^* Employing this method, in this work, 
we present the quantum spectra of the nuclear tetrahedral rotor, which are obtained 
microscopically for the first time, and discuss briefly its characteristic properties. 
More detailed investigations will be reported elsewhere. 
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§2. Efficient Method for Projection and GCM 

Since the detailed information is published in Ref. |7j), here we briefly discuss the 
essential points of our method. The projector is of the form P a = J g a (x)D(x)dx, 
where D(x) is a unitary transformation of symmetry operations, e.g., the rotation, 
and the main task is to calculate a quantity (^>\OD{x) \<P') for an arbitrary operator O 
and HFB type states \^>) and \<P') at the mesh points of integration over the parameter 
space (x). The general HFB state is defined by the quasiparticle operators, 

/9 fc |*> = (fc = l,2,...,M), $l = Y,[ U "°$ + V «< £ i]> ( 2-1 ) 

l 

where (cj,c>) are the creation and annihilation operators of the states within the 
original basis, for which the unitary transformation D is represented by a matrix 
D = (A'i), 

Dc\b^ = Y J Dnc\, (2-2) 

V 

Then, the quantity (<P\OD(x)\<I>') can be expressed in terms of the matrix elements 
of O, the transformation matrix D, (U, V) amplitudes in Eq. (|2-1|) and the similar 
amplitudes (U',V) for |^'). All the matrices have sizes M x M and each matrix 
manipulation requires 0(M 3 ) operations. This is the difficulty for the projection 
from the HFB state with a large model space; for example, if we take N™ x = 20 
spherical harmonic oscillator shells, the dimension becomes M = 3, 542. If the 
pairing correlation is neglected and the calculation is restricted to the HF state, it 
is shown that the matrix operations reduces to 0(MN 2 ) for one-body operators, 
where N is the number of particles and is at most 200 or so. 

The method we have taken to circumvent this difficulty can be summarized by 
the following two points: 

• Basis truncation in terms of the canonical basis; with 

tfc = Yl Wlkd l l # > = hk'vl (v 2 k in descending order), (2-3) 

the P space is defined as {k = 1, 2, L p (e); v\ > e} with a small number e. 

• Full use of the Thouless form with respect to a Slater determinant; 

i#> = n«p(2>«'"!4W i^rHio), 4 = { S 

V />Z' ' k=l k 

(2-4) 

In order to avoid the sign problem of norm overlap for general HFB statesp* espe- 
cially those without the time-reversal invariance, one has to calculate the pfaffian 
with Thouless amplitudes. Then it is easy to show that the matrix operations for 
one-body operators can be reduced to 0{ML^{e)) by the canonical basis truncation 
above. However, the Thouless amplitude with respect to nucleon vacuum diverges 
when the pairing correlation is vanishing or the HFB state is orthogonal to the nu- 
cleon vacuum, e.g., a quasiparticle excited state. The Thouless form with respect to 
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a Slater determinant in Eq. (|2-4p solves this difficulty because then the amplitude 
Z ~ Uk/ V k for k < N and Z ~ Vk/uj, for A; > N in the diagonal representation, 
and the truncation scheme can be utilized for any HFB type mean-field states. In 
addition, the calculation can be further reduced because of the simple treatment of 



the core space (deep hole states) defined as {k 

Examples of the occupation and 
empty probabilities, v\ and u\, are 
shown in Fig. [TJ and the resultant two 
dimensions, L p (e) and L a (e), as well 
as the calculated rotational spectra are 
shown in Fig. [21 It can be seen that the 
truncation parameter e ~ 10~ 3 — 1CP 4 
is sufficient, and then the necessary P 
space dimension is L p « 150 (L Q 50), 
which is much smaller than the original 
basis size M ~ 3, 000. Thus the calcu- 
lational effort is dramatically reduced. 



l,2,...,L (e);u( = l-v 2 k >e}. 
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Fig. 1. Probabilities v\ and u\ as functions of 
the number k of the canonical basis. 
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Fig. 2. Dimensions of the P-space L p and the core space L as functions of the model space trun- 
cation parameter e (left), and the rotational spectra (J = 2, .., 8) as functions of the truncation 
parameter e for 164 Er (right). 

As for the choice of the Hamiltonian, it may be desirable to employ effective 
interactions like Skyrme and Gogny forces, but they have problems for the projection 
and the GCM calculations due to their density-dependence. Thus, for the illustration 
of our method, we use the following schematic multi-separable interaction: 



H = h 



E 

A,r=n,p 



9 a CT\ ■ G 



(2-5) 



where the single-particle part h is composed of the Woods-Saxon potential. The two- 
body interaction is constructed from the isoscalar particle-hole channel operators 



^ r=n p F£ and the pairing channel operators G\ u defined by 



(r) 



dr 



A/i I 



G 



Ro 



4tt 



2A + 1 



Y Xlx (e,<t>), (2-6) 



with VJ(r) and Kk~ being the central part of the Woods-Saxon potential and its radius 
parameters, respectively, and Rq = 1.2A 1 / 3 fm. The multipolarities A = 2,3,4 
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are included in the p-h channel and A = 0,2 in the pairing channel. As for the 
force strengths, we adopt the selfconsistent value for while the values of g$ are 
determined according to even-odd mass differences A T . The ratio g^/do = 13.6 is 
chosen to reproduce typical rotational bands in a rare-earth and an actinide nucleus, 
see Ref. |7j) for details. It should be noted that this kind of schematic Hamiltonian 
cannot be used to obtain the ground state energy and deformation; it is only used 
to calculate the collective excitation spectra based on a given mean-field state \&). 



§3. Spectra of Nuclear Tetrahedral Rotor 



The tetrahedra deformation can be conveniently described by the 032 deforma- 
tion without other terms in the usual nuclear surface parametrization, 



R{0, V ) = iW{<*}) 



O 



1^(9, y>) 



(3-1) 



Thus, it is clear that both the axial symmetry and parity are broken. Now let us 
investigate what kinds of quantum spectra with definite spin-parity are expected 
for the tetrahedrally deformed nuclear state. For this purpose, it is interesting to 
consider the region near one of double tetrahedral-closed shell nuclei, 110 Zr (Z = 40 
and N = 70), because experiments for neutron-rich Zr isotopes have been done 
recently at RIKEN. 91 In this experiment a new isomer with life time of about 600 nsec 
has been found in a nucleus 108 Zr, and speculated as a possible low- lying tetrahedral 
state according to a theoretical prediction!^ 
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Fig. 3. Calculated spectra of the ground state band in 108 Zr (left), and of tetrahedral states 
assuming Q32 = 0.14 (right). Dotted curves denote the /(/ + 1) spectra going through 3J~. 
Although higher excited states may not have realistic meaning, they are included to show what 
kind of spin-parity appears. 

We apply our quantum number projection method to 108 < n 02 r nuclei. In order 
to obtain a reliable mean-field for the ground states, we have performed the paired 
Woods-Saxon Strutinsky calculation of Ref. [TT]) (with the Wyss-2 potential param- 
eter set). Since the pairing model space is cut off in these calculations, we do not 
include the cut-off factor, which was introduced in Ref. [7]) for the pairing operators 
G\ in Eq. (j2-5j) . in the present projection calculation. The pairing force strength g$ 
are fixed according to those results with the same ratio g^/do = 13-6 as previously. 
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The model space size N£ 



18 and the truncation parameter e = 10 4 are used. 



Because the cranking procedure is important to evaluate the moment of inertia^ we 
use a small cranking frequency fuv rot = 0.01 MeV about the a>axis in all calculations. 

We compare calculated ground state bands in 108 Zr with experimental datcP in 
Fig. [3l Since the total energies of prolate and oblate minima are almost the same, we 
show spectra on both of them. The oblate spectra better agrees with data, though 
the prolate ground state energy is slightly lower. In any case, this confirms that our 
Hamiltonian is a reasonable one. The calculated spectra of tetrahedral states with 
assumed deformation^ 1 (X32 = 0.14 are also shown in Fig.[3j Although the spectra are 
not so rotational-like, their spin-parity appears to be characteristic to the simplest 
singlet representation of the tetrahedral rotor 1 ^* corresponding to the closed shell, 
i.e., + , 3~, 4 + , 6 ± (doublet), etc., even for a neutron non-closed nucleus 108 Zr. This 
is because the totally symmetric state is realized owing to the pairing correlation. 
If the pairing is switched off, two neutrons partially fill the last fourfold generate 
orbit and the resultant Slater determinant is not tetrahedrally symmetric. Then, 
the spectra are more complex combinations of various irreducible representations. 
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Fig. 4. Calculated spectra of tetrahedral states in 110 Zr with various values of (J32. These results 
correspond to projection from a single HFB wave function; the high lying bands serve merely 
the illustration of the symmetry properties and are not meant to be compared with experiment. 

In Fig. H]the calculated tetrahedral spectra of the double closed nucleus 110 Zr, 
for which the mean-field state is tetrahedrally symmetric and gives the same char- 
acteristic spectra even in the unpaired case, are depicted for various values of 032- 
It is clearly seen that the spectra gradually change to rotational ones at consider- 
able deformation 0:32 = 0.25, while they are more like vibrational or transitional in 
0.1 < 032 < 0.20; the 3j~ energy quickly decreases as Q32 increases. This kind of 
phase transitions is well-known for the quadrupole deformation. We plot the moment 
of inertia estimated by 6/E{3^) in Fig. [SJ where the results with switching off the 
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Fig. 5. Calculated moment of inertia deduced from the 3j~ energy (left) and the B(E3 : 3]" Of) 
in Weisskopf unit (right) as functions of 032 for 110 Zr. The rigid-body inertia is 35.1 [ft 2 /MeV]. 

pairing correlation are also included. Note that the effect of the pairing correlation 
is not so strong. This may be partly because this nucleus is double tetrahedral- 
closed and the pairing gaps of neutron and proton are not so large. Compared to 
the rigid body value, ^AMRq = 35.1 h 2 /MeV, the moment of inertia is rather small 
even at 032 ~ 0.4 without pairing correlation. The calculated B(E3 : 3± — > Of) 
value is also shown in Fig. which well corresponds to the rotor model estimate, 
(2-5 K o)\{&\r 3 Y 3K \<P)\ 2 (3K3- K\00) 2 , with K = 2. Again, they become very large, 
more than 100 Weisskopf unit, for 032 > 0.3. 

In summary, we have calculated the spectra of the nuclear tetrahedral rotor 
for the first time by employing the quantum number projection method. By in- 
creasing the tetrahedral deformation, the spectra gradually change from vibrational 
to rotational just like in the spherical-to-deformed phase transition in the case of 
quadrupole deformation. It is quite interesting to measure the characteristic spec- 
tra of the tetrahedral rotor, which is a clear indication of the exotic tetrahedral 
deformation in atomic nuclei. 
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